Photon production from the vacuum close to the super-radiant transition: 
When Casimir meets Kibble-Zurek 
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The dynamical Casimir effect (DCE) predicts the generation of photons from the vacuum due to the para- 
metric amplification of the quantum fluctuations of an electromagnetic field. The verification of such effect 
is still elusive in optical systems due to the very demanding requirements of its experimental implementation. 
This typically requires very fast changes of the boundary conditions of the problem. We show that an ensemble 
of two-level atoms collectively coupled to the electromagnetic field of a cavity (thus embodying the quantum 
Dicke model), driven at low frequencies and close to a quantum phase transition, stimulates the production of 
photons from the vacuum. This paves the way to an effective simulation of the DCE through a mechanism 
that has recently found experimental demonstration. The spectral properties of the emitted radiation reflect the 
critical nature of the system and allow us to link the detection of DCE to the Kibble-Zurek mechanism for the 
production of defects when crossing a continuous phase transition. 

PACS numbers: 42.50.Pq, 32.80.Qk, 64.60.Ht 



When N two-level atoms interact collectively with a single 
mode of the electromagnetic field inside a cavity, thus realiz- 
ing the so-called Dicke model 1 1], a critical value of the atom- 
photon coupling Qc exists at which the system undergoes a 
quantum phase transition, generally referred to as the super- 
radiant transition. Below such coupling, the atoms are in their 
ground state and the cavity field is unpopulated. Conversely, 
above gc there is a spontaneous symmetry breaking and the 
photon field gets populated through a mechanism producing 
a displaced coherent state ^. The experimental demonstra- 
tion of the super-radiant transition in the Dicke model has re- 
mained outstanding until recently, when a key result has been 
achieved in a set-up involving intra-cavity Bose-Einstein con- 
densates [3]. A super-radiant transition has been enforced 
by exploiting the spatial self-organization of the atoms in an 
intra-cavity condensate coupled to the cavity field and sub- 
jected to an optical-lattice potential. 

Here we investigate the relation between equilibrium and 
dynamical properties of a Dicke system brought close to a 
quantum phase transition. We prove that, at the super-radiant 
transition, a DCE-like mechanism [4, 5] arises from the use of 
a time-dependent driving and results in a flux of photons gen- 
erated from the vacuum fluctuations. DCE has been predicted 
to occur in QED settings involving a cavity with oscillating 
end mirror [6] . This scheme, however, appears to be techno- 
logically demanding given the prohibitively large frequency 
at which the mirrors should vibrate to produce a measurable 
flux ofphotons. Notwithstanding some interesting propos- 
als IlZI-IsI] having the potential to ease the requirements for its 
observability, an experimental demonstration of DCE is still 
elusive in the optical domain. Recently, a DCE-like mech- 
anism has been observed in an experiment performed using 
microwaves iIiqIi . 



Our proposal pursues a different direction: we observe that, 
on approaching the Dicke super-radiant phase transition, the 
frequencies at which the DCE-like effect becomes observable 
are lowered, thus narrowing the gap separating the experimen- 
tal state-of-the-art from the observation of the effect. More- 
over, we unveil an intriguing connection between the occur- 
rence of DCE through the mechanism we propose and the 
Kibble-Zurek mechanism (KZM) [11, 12]. The latter pre- 
dicts the formation of defects in a quantum many-body system 
dragged through a critical point 11 13141511 and is due to the in- 
ability of the system to remain in its ground state. The produc- 
tion of defects occurs regardless of how slowly the dragging 
is performed and the mechanism has been shown to be related 
to adiabatic quantum computation [16] and quantum anneal- 
ing [17]. We are thus able to bridge two fundamental phe- 
nomena in out-of-equilibrium quantum systems with the goal 
of simplifying their observation. The recent demonstration 
of the Dicke super-radiant transition L3J, which is the build- 
ing block of our proposal, marks a promising starting point 
towards an experimental investigation along the lines of our 
work. 

We consider a cold system of N two-level atoms, collec- 
tively interacting with the field of a single-mode cavity whose 
annihilation (creation) operator is d(d^). Each two-level atom 
is modelled as a pseudo-spin whose Pauli spin matrices are 
{(j^, (J* } (i=l, .., N). We realistically assume a small-sized 
atomic sample, neglect the variations of the cavity field at its 
location and take the coupling strength as uniform. Then, we 
introduce the total angular momentum J of the atomic sample 
with components J± = ^^ a]_ and Jz = J]* ^l ^^^ ^^^~ 
sider time-dependent single-spin energy splittings uJi){t). The 
Hamiltonian of the system in the dipole approximation thus 



0.006 


2eo 

1 


0.005 


0.004 






0.003 






0.002 






0.001 


Jv 



{a.^ a) 




FIG. 1 . Radiation flux outside the cavity, (a) Flux of photons outside 
the cavity against 77 for g/uJa = 0.9, Qc/ooa = 0.45, ^/uJa = 0.005, 
and \/uJa = 0.005. For these parameters, eo/cJa ~ 0.315. (b) 
Flux of photons outside the cavity against 77 and g for ujh/ooa = 1, 
7/cJa = 0.005 and X/ua = 0.005. 



reads (we set /i = 1 throughout the paper) 

^0 = ^aCi^ci + uJh{t)Jz + -^{o) + a){Jj, + J_), (1) 



2j' 

where cja is the frequency of the cavity and cf is the atom- 
field coupling constant. Analogously to Refs. 1 81191], we take 
ujh{i) = cjo + Asin(7^t). The parameter j is the cooperation 
number in the Dicke theory, that is an eigenvalue of J^ which, 
together with the eigenstates of J^, is used to build the Dicke 
states. The ensemble of N two-level atoms is then described 
as a pseudo-spin of size j = N/2. In this case, photons are 
generated at r^res = 2co'a. In Ref. |18] a model similar to ours 
but based on a semi-classical approach has been addressed 
to relate DCE-like effects to Dicke super-radiance. Here we 
perform a full quantum treatment of both the atom-light in- 
teraction and the effects on the photon statistics induced by 
the driving of the atomic subsystem. Moreover, as discussed 
in the second part of the paper, we will unveil the connection 
between the DCE-like effects and the KZ mechanism. 

Eq. (O strongly resembles the many-body Landau-Zener 
problem studied in Refs. [19, 2QJ. A crucial difference be- 
tween the two cases is the presence of the counter-rotating 
terms in Eq. Q. These lead both to the super-radiant transi- 
tion and the production of photons. For N :^ l^thc Holstein- 
Primakoff representation of the angular momentum 112 IH can 
be used to approximate the atomic cloud to a non-linear har- 
monic oscillator. However, for a large atomic sample (as in 
the case here), we can take J+ ^ V^jb^ with [6, 6^] = 1 and 
retain the harmonic approximation. Eq. ^ thus reduces to 

H = ujad^d + ujb{t)b^b + g{d^ + d){b^ + 6), (2) 

which is easily diagonalised: the normal frequencies e±{t) 
and modes q± {t) are discussed in Appendix A (their form is 
not relevant for our purposes). The normal-mode description 
of Eq. (O simplifies the analysis of the critical properties of 
the time-modulated Dicke model [22], which shows the exis- 
tence of a critical value gc (t) = ^/iJaUO^(f) / 2 at which a phase 
transition occurs. In the phase corresponding to g < g^ which 
we dub normal phase, the number of photons in the cavity 
mode is very small and (a) = (b) =0. Upon approaching 



gc{t), the number of photons increases and the system reaches 
the so-called super-radiant phase at g>gc{t), where the cav- 
ity mode is macroscopically populated even in the ground 
state of the system and (k) ^ {k=a^b) due to spontaneous 
symmetry-breaking: while at small couplings the total number 
of excitations {d^d + b^b) is conserved, at the phase transition 
such a symmetry is broken, resulting in the population of the 
cavity field. Here, we will only refer to the normal phase. 

Let us first address the case of a lossless evolution. For pa- 
rameters of the system such that e+(t) :^ {e_ (t), r]}, the non- 
critical mode ^+ will not contribute to the photon production 
and we can consider only q- . To simplify the notation, we 
drop the index and set e_(t) = e{t) from now on. Moreover, 
taking A small, the system can be treated as a harmonic os- 
cillator with frequency eo perturbed by a weak driving at the 
modulation frequency 77. In this framework, excitations can be 
created only at the resonance condition r] = keo (k e 1^). If 
the system is initially in its ground state, which corresponds to 
the vacuum of the effective harmonic oscillator, such dynam- 
ics lead to the generation of photons from the vacuum inside 
the cavity, along the lines of DCE, already at r] = 2eo. The 
number of photons generated increases as ^ -^ gc. This is 
very important: as eo ^ when the system approaches its 
critical point, the r] needed to observe the DCE-like effect is 
lowered, thus bringing its verification closer to experimental 
feasibility. 

The assumption of unitarity is not realistic and we now in- 
clude the leakage of photons from the cavity, a mechanism 
that is conveniently tackled by means of the input-output for- 
malism for optical cavities ll23ll and by modelling the bath 
to which the cavity field is coupled as a distribution of har- 
monic oscillators with associated operators {a^^al), having 
frequency u and interacting with the cavity field with strength 
kjy according to the Hamiltonian V=i J^ kj,(Qf^a^ — dal)du. 
The method used to tackle such open dynamics is fully de- 
scribed in Appendix B. Here we only state that, by assum- 
ing a flat density of states in the bath p(i/) vanishing for 
uKO and kjy=k (Vi^>0), manageable expressions for the ef- 
fective dissipation rates affecting the evolution of the sys- 
tem at hand are found. Moreover, we can determine the sta- 
tionary mean number of photons inside the cavity (a^a) = 
liint^Qo{a\t)a{t)) and, from this, the out-coming photons 
flux (oftd) = /Q'^dcj«^*^<^*). Here, <^^ (d^^^t) is the 
field annihilation (creation) operator of the output mode at fre- 
quency z/ 1I24I] . 

We can now discuss the qualitative features of the mech- 
anism here achieved. First, we find no out-coming photon 
for an unmodulated driving field. On the contrary, any time- 
dependent modulation generates a constant flux of photons. 
This is clearly understood by taking the case of a small- 
amplitude modulation and stopping at the first order in A ll25ll . 
In Fig. [U we show the out-coming photon-flux for a damp- 
ing rate 7 (positive and constant for a;>0) against r] and g. 
In Fig. [T](a) a resonance peak is clearly visible at 7^ ^ 0.63 
when g = 0.45, corresponding to ^ = 0.9 gc and eo ~ 0.315, 
which confirms that a resonance peak is achieved at 7/ ^ 2eo. 
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FIG. 2. Spectral density of the output photons. Taking c^a = c^b = 1, 
A = 0.005, 7 = 0.005, g = 0.9 and Qc = 0.45, we find eo = 0.315. 
We have taken r]/2eo = 1 (corresponding to resonance conditions, 
main panel), r]/2eo = 0.7 (upper inset), r]/2eo = 1.3 (lower inset). 



Such a prediction is strengthened by Fig.[T](b), where the flux 
of out-coming photons is plotted against g and r] and shows 
resonances at t^ ^ 2eo, regardless of the value of eo. 

We complete our study by considering the spectral den- 
sity S{uj) = (a°^*^(a;)d°^*(cL;)) of the output field, which is 
plotted in Fig. [2] for a weak modulation driving at variable r]. 
When 7/=7^res=2eo, the spectrum reveals a single sharp peak 
at cj ^ eo [cf. Fig. [21 main panel]. In the non resonant regime, 
the emission Situj ^ cq is drastically reduced and sidebands of 
enhanced emission atuj ^ r]±eo and u ^ r]-\- 2eo appear (cf. 
insets of Fig. O. 

A series of remarks are due. First, we reassure on the va- 
lidity of the Holstein-Primakoff approximation. In all our cal- 
culations of the out-coming photon-flux, the number of exci- 
tations in the atomic system is much smaller than the actual 
number of atoms in the sample, thus ensuring that our frame- 
work holds. Second, due to the small frequency of the pho- 
tons generated when g approaches the critical value (the emis- 
sion frequency is cj ^ eo), thermal noise in the output signal, 
which has not been included in our study, may be significant. 
However, the generated photons can be detected also in the 
presence of strong background noise simply by using a cav- 
ity with semi-transparent mirrors. In this scheme, the cavity 
mode is coupled with two thermal baths and photons are al- 
lowed to enter/abandon the cavity from both sides. A the two 
baths are uncorrected, the noise can be virtually eliminated 
by measuring the correlations between the output modes. 

Finally, we address the crucial connection between our 
DCE-like mechanism and the Kibble-Zurek one lllll Il2n . On 
approaching the critical point of the model in Eq. ([3]), regard- 
less of the value of r], there will always be a regime where 
the perturbation is non-adiabatic and photons are produced. 
A first estimate of the unavoidable departure from adiabatic- 
ity, with a consequent photon-flux, is obtained by calculating 
the probability of the system to go into an excited state. For 
simplicity, we consider one period in the absence of damp- 
ing. The probability of leaving the ground state at the final 




FIG. 3. (a) Schematic representation of the four freeze-out points in 
the trigonometric circle, (b) Probability of leaving the ground state 
against 77/ eo for g = 0A9/(jJa and various values of A. 



time tf (ti being the initial time) is P=l-\{'^{tf)\(po{tf))\'^ 
with \(pn{t)) the instantaneous eigenstates of the harmonic os- 
cillator and |^(t/)) the final state of the system. The KZM 
relies on the assumption that the state of a system brought 
close enough to the critical point freezes when the system is 
not able to adiabatically follow the changes in the control pa- 
rameter. For the driving here at hand, the freeze-out times 
is found by solving the equation T{t)/T{t) = r(t), where 
T{t) = gc{t)/g — 1 plays the role of the relative temper- 
ature of the system and r = ro/e{t) is its relaxation time 
(to = l/uj) 1 26]. For a sinusoidal modulation of T{t) and 
if the oscillating terms brings the system sufficiently close 
to the critical point, one finds four solutions, each embody- 
ing a freeze-out time. Fig. [3] (a) shows their representa- 
tions in the unit circle. As the system is initialized in its 
ground state, i.e. |^(ti)) = \(po{ti))^ the adiabatic condition 
T{t)/T{t) > r is satisfied until t = £1, where £1 is the freeze- 
out time at which the system enters the so-called impulsive 
regime. During this period, the state of the system is frozen 
until t = £2, when the adiabatic condition is fulfilled again 
and the state of the system becomes |^(£2)) = |^o(£i)) = 

Zln ^^,o(^^2,£l)|^n (£2)) with Cn,m(^,tO = {^n{t)\(Pm{t^)) . 

The same argument applies to the second part of the cy- 
cle, where the system evolves adiabatically for t G [£2, £3] 
and is frozen for t G [£3, £4]. The state at £4 is then 

l^(^^4)) = E/c,nCfc,n(£4,£3)Cn,0(£2,£l)e"'^H^/c(£4)) with 

On = Jj:^ dtEn{t). Finally, the last part of the evolution 
(t e [£4, t/]) will not affect the probability P, which is thus 
P = 1 — |(^(£4)|(/:)o(£4))P and whose behavior against r] is 
shown in Fig.[3](b) for different values of A. Clearly, the closer 
the system to the quantum phase transition, the more it is sus- 
ceptible to a low-frequency driving. A more detailed analysis 
requires the study of the transient dynamics. The scheme of 
Fig. [3] (a) is still valid, the probablity of excitations being cal- 
culated by composing four different dissipative maps in the 
same spirit of Ref. 112711 . We only expect quantitative changes. 
To corroborate the connection between DCE and KZM, we 
have further analyzed the photon production in the adiabatic 
and non-adiabatic regimes [cf. Fig. IH. For r]>emin (being 
emin the minimum value of e{t) over a cycle), the dynamics 
is non-adiabatic and photons can be created. Close to criti- 
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FIG. 4. Output photon-flux as a function of 77 for different values of 
g. The transition between adiabatic and non-adiabatic regime (sharp 
step) is located at the minimum of the gap and is shifted to lower 
frequency when the coupling gets closer to the critical coupling. At 
the critical point the dynamics is purely non-adiabatic. 



cality, the minimum of the gap vanishes, the system is always 
in the non-adiabatic regime and the photon-flux increases lin- 
early with ?7 until the maximum value at resonance is reached. 
Far from the transition, the photon production decreases from 
the resonance with a Lorentzian behavior: when 7/<emin, the 
photon-flux is sharply reduced and a linear behavior is recov- 
ered but with a much smaller value. This abrupt transition be- 
tween the adiabatic and non-adiabatic regimes demonstrates 
that the breakdown of adiabaticity due to critical slowing- 
down is at the origin of photon creation in the DCE, a situation 
totally analogous to what is described by the KZM. 

We have proposed a scheme to achieve DCE-like effects 
by dragging a driven quantum Dicke model across its critical 
point. By linearizing the model, we have related the mech- 
anism of photon generation from the vacuum to the proper- 
ties of the eigen-modes of the system, thus providing a clear 
picture of the DCE-like effect arising from a Dicke quan- 
tum phase transition guided by a time-modulated driving. As 
the frequency of the driving is quenched at the Dicke criti- 
cal point, the observation of a sizable flux of generated pho- 
tons becomes less demanding for an optical-domain verifica- 
tion. We have also connected the photon-generation process to 
the KZ predictions for defect nucleation at a quantum critical 
point, thus pointing out the intimate connection among three 
fundamental mechanisms in quantum many-body physics. 

We acknowledge financial support by the National Re- 
search Foundation and Ministry of Education in Singa- 
pore, the EU FP7 Programme under Grant Agreement No. 
234970-NANOCTM and No. 248629-SOLID, UK EPSRC 
(EP/G004759/1), and EUROTECH. 



[6] V. V. Dodonov, Phys. Scr. 82, 038105 (2010). 

[7] J. R. Johansson, G. Johansson, C. M. Wilson, and F. Nori, Phys. 

Rev. Lett.103, 147003 (2009). 
[8] A. V. Dodonov, etaL, arXiv:0806.4035Y3 (2008). 
[9] S. De Liberato, D. Gerace, I. Carusotto, and C. Ciuti, Phys. Rev. 
A80, 053810 (2009). 

[10] C. M. Wilson, G. Johansson, A. Pourkabirian, M. Simoen, J. R. 
Johansson, T. Duty, F. Nori, and P. Delsing, Nature (London) 
479,376(2011). 

[11] T. W B. Kibble, J. Phys. A 9, 1387 (1976). 

[12] W. H. Zurek, Nature (London) 317, 505 (1985). 

[13] J. Dziarmaga, Adv. Phys. 59, 1063 (2010). 

[14] W H. Zurek, U. Dorner, and P Zoller, Phys. Rev. Lett.95, 
105701 (2005). 

[15] A. Polkovnikov, Phys. Rev. B72, 161201(R) (2005). 

[16] E. Farhi, et ai, Science 292, 472 (2001). 

[17] G. E. Santoro, etaL, Science 295, 2427 (2002). 

[18] L L. Kirilyuk, Opt. Spectrosc. 92, 719 (2002). 

[19] A. Altland, V. Gurarie, T. Kriecherbauer, and A. Polkovnikov, 
Phys. Rev. A79, 042703 (2009). 

[20] A. AMand, V. Gurarie, Phys. Rev. Lett.lOO, 063602 (2008). 

[21] T. Holstein, and H. Primakoff, Phys. Rev. 58, 1098 (1940). 

[22] C. Emary, and T. Brandes, Phys. Rev. E67, 066203 (2003). 

[23] D. F. Walls, and G. J. Milburn, Quantum Optics (Springer, Hei- 
delberg, 2008). 

[24] C. Ciuti, and L Carusotto, Phys. Rev. A74, 033811 (2006). 

[25] A. Mari, and J. Eisert, Phys. Rev. Lett.103, 213603 (2009). 

[26] Notice that T{t) is defined as T(t) = 1 - g{t)/gc if the time- 
dependent parameter is g instead of cj. 

[27] A. V. Shytov, D. A. Ivanov, and M. V. Feigel'man, Eur. Phys. J. 
B 36, 263 (2003). 



APPENDIX A: FREQUENCIES OF THE NORMAL MODES 

In this Section we provide additional details on the achieve- 
ment of the two-oscillator model used in the main Letter and 
give the explicit form of the normal-mode frequencies. We 
start referring again to the Hamiltonian of the system in the 
dipole approximation reads (we set fi = 1) 



9 



Ho = ujad'^d + uJb{t)J^ + ^(^ + ^)(^+ + J-)^ (3) 



2j 



where, as discussed in the main Letter, Ua is the frequency of 
the cavity and g is the atom-field coupling constant. For N :^ 
1, the use of the Holstein-Primakoff picture is well motivated. 
The angular momentum describing the atomic cloud is thus 
reduced to a non-linear harmonic oscillator as 






(4) 
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Here, h and ¥ are the operators of the effective non-linear 
boson. As the number of atoms is large, j ^ 1 and a 
good approximation is given by taking J+ = Jl ?^ y/^o^ • 
The Hamiltonian thus reduces to Eq. (2) of the main Let- 
ter, which is straightforwardly diagonalized to give the eigen- 



frequencies 

' (5) 

These are associated with the eigen-modes q± {t) whose ex- 
pHcit form is not relevant for our analysis. 



APPENDIX B: LOSS-AFFECTED DYNAMICS, 
OUT-COMING PHOTON FLUX AND SPECTRUM 



In the losses-affected scenario, the dynamics of the cavity 
mode is conveniently tackled by means of the input-output 
formalism for optical cavities ll23ll and by modelling the bath 
as a continuous distribution of harmonic oscillators with fre- 
quencies v, formally described by the bosonic operators d^^ 
and al and each coupled to the cavity field with a strength 
kjy. As in Eq. (3) of the main Letter, the cavity field- 
bath coupling Hamiltonian is V= J^ di'kiy{aiya^ — aaj,). 
The dynamics of the system is analysed by studying the 
Langevin equations [Is], which can be cast in a compact ma- 
trix form. Taking into account that the effective atomic mode 
b does not experience dissipation for a fully condensed gas 
without atom losses, we define the bosonic operators vector 
u{t) = (a(t), 6(t), d^t)^ ^U^))^ ^^^ the Langevin-force vec- 
tor F{t) = {fit), 0, p(t), 0)^, such that the equations of mo- 
tion for the atomic and cavity-field modes read 



ii{t) = -iM{t) ■ u{t) - / dt'T{t - t')u{t') + Fit). (6) 

Here M(t) is a time-dependent 4x4 matrix 
that takes into account the unitary evolution and 
Vit-t')=&idig[-f{t-t'), 0, ^{t-t'), 0] is the dissipation kernel 
whose elements j{t — t') are decay rates. The explicit form of 
such rates will be discussed later on. In light of our choice for 
the modulation of the atomic energy splitting, the elements 
of M{t) are oscillating functions with period T = 27r/r], so 



that we can take M{t) = Mq 
time-independent term 



Mi(e*^^-e-*^^) with Mq the 



f^a 



Mn = 



9 







g ^0 g 

-g -UJa 

\-g -g 



3 \ 



-g 



(7) 



and Ml = diag[0, A, 0, —A] that takes into account the mod- 
ulation. Notice that, due to the inclusion of the counter- 
rotating terms, the decay rates mVit — t') in Eq. (O are time- 
dependent. Solving Eq. (O is made difficult, in general, by the 
convolution integral describing damping. The approach is sig- 
nificantly simplified by moving to the frequency domain. The 
Langevin-force operators f{uj) in the frequency domain are 
linked to the input noise operators as /(cj) = 2Tik^piuo)a^^, 
where p{uj) is the density of states of the bath. The decay 
rates 7(co') G C are thus found to have 5R[7(cl;)] = 7r|/Cc^ pp(cj) 



and ^[7(cj)] = -^V j^^duj'^[^{uj')]/{uj' - cj), where V 
denotes the principal value of the integral f24]. While the 
imaginary part of 7(0;) is just a frequency shift, 3?[7(cj)] is 
responsible for the frequency-dependent damping of the cav- 
ity mode. When the counter-rotating terms are taken into ac- 
count, it becomes crucial to consider that the density of pho- 
tonic state in the bath p(cj) is zero for negative frequencies. 
Indeed, by modeling the bath as a collection of harmonic os- 
cillators, only positive frequencies have a physical meaning. 
Within this assumption, it follows immediately from the def- 
initions given above that 7(0;) = and f{uj) = for cj < 0. 
We will also suppose that the damping rate is constant for pos- 
itive frequencies. This is equivalent to assuming that k^ = k 
and, for a; > 0, p{uj) = 1 Within this assumption, we de- 
fine 70 = 7r|A:p and we can write Re[j{uj > 0)] = 70 and 
Re[7(a; < 0)] = 0. 

In the following the index m = 0, ±1 indicates the number 
of sidebands in the drive and Qij (uj) are the matrix elements 
of ^(cj) = [iM{uj)]~^. We first consider the mean number 
of photons inside the cavity at the stationary state (d^d) = 
liint^cx){d^ {t)d{t)) . After Fourier transforming d{t) and in 
the limit t^oo we have (a''"a) = (27r)~^ f^ duj{d^{uj)d{uj)). 
Using the steady state solution of the cavity operator and 
the expression for f{uj) given in the last paragraph, fiuj) = 
27rk^p{uj)a]^, the mean number of photons inside the cavity 
at the steady state is given by 
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dc^|S4m+l,4(m+j)+3(-^)l • (8) 
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Here we have assumed that the input field is in the vacuum 
state, so the operators d^J^ fulfill the condition (ofjj^ajy ) = 
S{uj — uj'). A factor of p{—uj — jr]) is responsible for the shift 
in the integration limits. 

In order to obtain the output operator of the cavity a°^*, 
we substitute the solution for the cavity operator d{uj) into the 
input-output relation [24]. Doing so, the output operator reads 



as 



j=-m 
m 

+ 270 ^ S4m+l,4(m+j)+3(^)p(-^-J^)Q^-u;-jr;- (9) 



j = -m 



The expression for the photonic flux outside the cavity is 
(oftd) = /Q^dcj(d^^*^a^^*). Since p{uj) = Oforcj < 
the only non- vanishing terms in the summation are the ones 
with j < 0. Notice that the negative frequencies have been 
taken out from the integration since only positive frequencies 
are physically allowed in the output field. Having the matrix 
elements of Q the dimension of time, the quantity (d^d) has 
the correct dimension of 1/t for a flux of photons. 

Finally, in order to calculate the spectral density S{uj) of 
out-coming photons, we first evaluate the analogous quantity 
V{uj) = {a\uj)a{uj)) inside the cavity. From this we find 



S{u) using the input-output relations, with the result S{u) = {jo / 7r)V (uj) . 



